A new insight into the observation of spectroscopic strength reduction in atomic 
nuclei: implication for the physical meaning of spectroscopic factors 
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Experimental studies of one nucleon knockout from magic nuclei suggest that their nucleon orbits 
are not fully occupied. This conflicts a commonly accepted view of the shell closure associated with 
such nuclei. The conflict can be reconciled if the overlap between initial and final nuclear states in a 
knockout reaction are calculated by a non-standard method. The method employs an inhomogeneous 
equation based on correlation-dependent effective nucleon-nucleon (NN) interactions and allows the 
simplest wave functions, in which all nucleons occupy only the lowest nuclear orbits, to be used. The 
method also reproduces the recently established relation between reduction of spectroscopic strength, 
observed in knockout reactions on other nuclei, and nucleon binding energies. The implication of 
the inhomogeneous equation method for the physical meaning of spectroscopic factors is discussed. 
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The concept of magic numbers of neutrons and protons 
making up a nucleus is fundamental to our understand- 
ing of a wide range of phenomena from the properties 
and binding energies of nuclei themselves to the relative 
abundance of elements in the universe [l[ . The observed 
magic numbers are usually explained by a model (the 
Shell Model) in which independent nucleons fill single 
particle energy levels in a mean field according to the 
Pauli exclusion principle. Such a picture is similar to elec- 
tronic structure of atoms responsible for organising the 
chemical elements into the Periodic Table. It has been 
found, however, that the cross sections of the (e, e'p) re- 
actions on the closed shell nuclei 16 , 40,48 Ca, 208 Pb are 
50-60 % smaller than those expected from the indepen- 
dent particle model Q. Direct reaction theories of the 
(e, e'p) reaction predict that its cross section depends on 
the Spectroscopic Factor (SF) which is a measure of the 
occupancy of single proton levels in the target nucleus. 
The observed reduction of SFs appears to contradict the 
traditional view of 16 , 40,48 Ca, 208 Pb as doubly magic 
nuclei. 

Away from the closed shells, nuclei such as 7 Li, 12 C, 
30 Si, 31 P, 51 V and 90 Zr also reveal a similar reduction 
of SFs as compared to prediction of the OUlo shell model 
Q. The SF reduction is also found for other nuclei in 
a recent analysis of the (d, p) and (p, d) reactions 0] , in 
which the bound state wave functions of the transferred 
neutron are fixed by modern Hartree-Fock calculations 
and have shapes similar to those derived from (e, e'p). 

Recently, SF studies with radioactive beams have re- 
vealed a new phenomenon. It turned out that reduction 
of experimental SFs S exp , detemined as ratios of the mea- 
sured to theoretical cross sections, from the theoretical 
values Sth, obtained in the shell model, depends on the 
separation energy of the removed nucleon and on the nu- 
cleon type. It has been also discovered that the SF re- 
duction factor R s = S exp /Sth is concentrated around a 
straight line when plotted as a function of the difference 



between proton (S p ) and neutron (S n ) separation ener- 
gies, AS*, taken as S p — S n and S n — S p for proton and 
neutron knockout, respectively 

It is known that Sth agrees better with S exp if the 
model space, in which Sth is calculated, is increased, or 
in other words, if particle-hole excitations are allowed. 
Thus, a six-shell treatment of 16 shows that the per- 
centage of the OHuj component in it is ^48-60% @ and 
that the 16 SF changes from the Ohio value of 2 to 1.65 
when the model space increases to Ahw jjjj]. However, it 
is still higher than the (e, e'p) value of 1.27(13) Q sug- 
gesting that more major shells should be added to the 
model space, which contradicts the view of 16 as a dou- 
ble magic nucleus. The contributions from missing model 
spaces can be recovered by using correlated wave func- 
tions in ab-initio approaches. Indeed, the 7 Li proton SF 
calculated in the Variational Monte Carlo (VMC) agrees 
very well with S exp from (e, e'p) [8j. However, for ' Li, 
8 B and 9 C the SF reduction obtained by VMC calcula- 
tions is not sufficient (see Table I). Also, the ab-initio 
calculations are feasible only for light nuclei while the SF 
reduction is observed for nuclei as heavy as 208 Pb. 

In this letter, I show that it is possible to reconcile the 
double magic nature of 16 O with the observed 60% reduc- 
tion of its spectroscopic strength and at the same time 
to explain the observed R S (AS) dependence if an alter- 
native method to calculating SFs is used. This method 
allows minimal shell model spaces to be used and ac- 
counts automatically for excluded orbits. It can be ap- 
plied to any nucleus and can be introduced into existing 
shell model codes including those used by the community 
of nuclear experimentalists studying one nucleon removal 
reactions. Below, I present this method, emphasize its 
importance for explaining the phenomenon of SF reduc- 
tion and present numerical results for A < 16 nuclei. 

The theoretical SF for one- nucleon removal, Sy, is de- 
fined in a model independent way as the norm of the ra- 
dial overlap function Iij (r) with orbital momentum I and 
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angular momentum j, calculated between the wave func- 
tions \&j B and ^j A of two neighbouring nuclei B = A— 1 
and A with the total spin J B and J a'- 
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All available shell model codes calculate S^- from izj(r) 
obtained by direct evaluation (DE) of Eq. ([I]) using some 
model wave functions in truncated model spaces. The 
input to these shell model calculations includes matrix 
elements of the effective nucleon-nucleon (NN) interac- 
tion fitted to a range of nuclear spectra. They carry no 
information about the radial shapes of Iij (r) , crucial for 
calculating one nucleon removal cross sections. In most 
applications, these shapes are found from the separation 
energy prescription, not related to the shell model NN 
matrix elements. 

An alternative method to calculate hj{r) is to solve 
the inhomogeneous equation (IE) 

{* B \T A -f B -E A + E b \^a) = (V B \V B - V A \V A ), (2) 

originally introduced by Pinkston and Satchler [9j . Here 
Tj and Vi are the kinetic and potential energy operators 
while Ei is the total energy of nucleus i. The r.h.s. of 
© is treated as known. Eq. @ generates 2y(r) which 
automatically have the correct asymptotic shape, a fea- 
ture crucial for transfer reactions but not so for binding 
energy calculations. Earlier explorations of this method, 
reviewed in [Iol | , were based on separating the mean field 
part out of Vi and keeping only the valence nucleon space. 
They gave little information of utility of the method and 
were abandoned before 1980s. Later, a different strat- 
egy, applied in [ll|, [T3| to calculating the source term 
(^b\Vb — Va\^a), resulted in SFs different from tradi- 
tional shell model values. Neither the Pinkston-Satchler 
approach nor that of Refs. 11, IjJ have been considered 
in the context of the SF reduction phenomenon as both 
were used at the time when the R S (AS) dependence was 
not known. Here, I prove the legitimacy of the method 
of Ref. [13] and show its relevance to the SF reduction. 
According to [12], the solution of Eq. ||2J| is 
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Z B Z N e 2 /r and V x = Y^<j In both Eqs. © and ©, 
r (r') is the distance between the centre-of-mass of B and 
the removed nucleon, and Yj is the spherical function. 
The advantage of ([3|) is that it guarantees the correct 
asymptotic form for lf E (r) when the experimental value 
of e is used, whatever tyj B and \P j A are. 

Eq. (J3]) was obtained assuming that &j B and ^ j A are 
exact solutions of the many-body Schrodinger equation 
and that V contains bare realistic NN interactions. In 
this case, I® E (r) and l[ E (r), and the corresponding SFs 
S^ E and Sf E , should be equal. However, usually &j B 
and \1/ j A are replaced by model wave functions in trun- 
cated spaces. This raises the question about what should 
be used for V. To answer it, I consider an exact nuclear 
wave function 'J constructed from an uncorrelated state 
<&, defined in some truncated model space, for example, 
using the Unitary Correlation Operator Method [Hi ]: 



|*) = C|*)=e*p{-i£«fc}|*). 



(5) 



i<3 



Here C is the unitary correlator designed to shift nucle- 
ons away from each other whenever their uncorrelated 
positions are within the repulsive NN core. <E> is found 
from an effective Hamiltonian that contain effective in- 
teractions V consisting of V = C^VC and the terms 
arising from the kinetic energy operator [l3j . If wave 
functions from Eq. ([5]) are used in Eq. ([3]), then 



(#b|V|*a> = ($b\C b (V a - V b )C a \$a) 
= ($b\VnC nb \$a) = ($ s |f cfr |^), 



(6) 
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where integration over r' is implied and G;(r, r') is the 
Green function for a bound nucleon in the field of a point 
charge Z B corresponding to the momentum in, 

Gi(r,r>) = -^ e ~ lii ^F l (iKr < )W_ tl)l+i (2Kr > ). 

ft K 2 

(4) 

Here K = (2/xe/Ti 2 ) 1 / 2 , e = Eb — Ea, (J- is reduced mass, 
rj = Z B Z^e 2 n/h 2 k, F is the regular Coulomb function 
and W is the Whittaker function. Also, V = V A — Vb — 



where Ca = C B C NB , C NB = exp{-i YU=i 9in} and 
V/v = Tld = \ v iA, assuming for simplicity that Coulomb 
interaction is absent. Eq. (J6j> tells us that the effec- 
tive interaction V cS that approximates V when mod- 
elling ll E (r) using uncorrelated model functions $ b and 
$,4, differs from the effective interaction V that gen- 
erates them. Moreover, $s and $a depend only on ma- 
trix elements (tp ai (ri)ip a2 (r 2 )\v cS (r 12 )\ip a3 {r 1 )ip ai (r 2 )) 
in a chosen truncated space, where i/j a (r) is a single- 
particle wave function in the state a. Hence, 
l( E {r) and Sl E depend on them as well. But in 
addition, they also depend on matrix elements of 
V eS , (G i (r,r 1 )/(r 1 r)^ 2 (r 2 )|C cff (r 12 )|^ 3 (r 1 )^ Q4 (r 2 )) 
(if centre-of-mass motion is neglected), that carry infor- 
mation about coupling to missing model spaces. This 
conclusion follows from the Green function expansion 
onto complete set {ip a (r)}, which includes states from 
both truncated and missing spaces. Thus, these matrix 
elements are not constrained by binding energy calcula- 
tions and must be constrained by some other means. 

A quantity that can serve as a reference to callibrate 
V cS is the Asympotic Normalization Coefficient (ANC). 
It determines the magnitude of the Iij(r) tail [14[, de- 
pends on the same operator V cS and can be determined 
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FIG. 1: The overlap functions rI IE (r) calculated for M3YE 
for (a) (d\p), (b) ( 3 H|d), (c) ( 4 He| 3 He> and (d) ( 16 0| 15 N) in 
comparison to (a) the realistic deuteron wave function ob- 
tained with AV18 potential in [is! ]. (b,c) ab-initio overlaps 
obtained with AV18+UR( or UIX) interaction in [l9l.l2C| and 
(d) the overlap function derived from (e,e'p) reaction 



TABLE I: S IE = S;? /2 + Sjg /2 calculated with M3YE and 
RM3YE in comparison to S DE , experimental values Ss XV 0, 
HH HI Hi HI and ab-initio VMC SFs S ab @, EJ HI, Hi ■ 
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M3YE 


RM3YE 


Sexp 


Sab 


3 H 


d 


1.5 


1.21 


1.33 




1.30 


3 He 


d 


1.5 


1.22 


1.35 




1.32 


4 He 


3 He 


2.0 


1.29 


1.42 




1.50 


7 Li 


6 He 


0.69 


0.28 


0.33 


0.42(4) 


0.42 


7 Li 


6 Li 


0.87 


0.44 


0.46 


0.74(11) 


0.68 


8 Li 


7 He 


1.02 


0.38 


0.44 


0.36(7) 


0.58 


8 Li 


7 Li 


1.14 


0.65 


0.77 




0.97 


8 B 


7 Be 


1.14 


0.78 


0.91 


0.89(7) 


0.97 


9 Li 


8 Li 


1.04 


0.60 


0.70 


0.59(15) 


1.14 


9 Be 


8 Li 


1.13 


0.45 


0.49 




0.73 


9 C 


8 B 


1.04 


0.71 


0.82 


0.77(6) 


1.14 


10 Be 


9 Li 


1.93 


0.81 


0.88 




1.04 


10 Be 


9 Be 


2.67 


1.48 


1.68 




1.93 


12 B 


n B 


0.99 


0.97 


0.84 


0.40(6) 




12 C 


U B 


2.85 


1.55 


1.76 


1.72(11) 




i3 G 


12 C 


0.63 


0.63 


0.51 


0.54(8) 




14 C 


13 c 


1.87 


1.82 


1.49 


1.07(22) 




14 N 


13 N 


0.72 


0.60 


0.53 


0.48(8) 




15 N 


14 N 


1.48 


1.31 


1.06 


0.93(15) 




16q 


15 N 


2.13 


1.57 


1.29 


1.27(13) 





from peripheral transfer experiments. In 15|, the ver- 
tex constants, related to the ANCs by a trivial relation 
14j , were studied for Op-shell nuclei in the oscillator 0Tlui 
shell model. It was found that reasonable agreement be- 
tween measured and calculated vertex constants can be 
achieved if a version of the M3Y potential, constructed 
in [l6| to fit the oscillator matrix elements derived from 
the NN scattering phase shifts, is used for V cS . Below, I 
use this interaction (labled as M3YE) to calculate S{ E '. 
I show that they are reduced with respect to SR and, 
at the same time, are closer to experimental SFs. 

First of all, I test the method as applied to the well 
understood A=2 system, for which l[ E (r) is the deuteron 
wave function and satisfies 



dr'r'Go(r,r')V cS (r')tp as (r'), (7) 



where ipo s is the 0s oscillator wave function. rlQ E (r), 
calculated with M3YE for V cS and with oscillator ra- 
dius r osc =1.51 fm, is close to the realistic deuteron wave 
function generated by the NN potential AV18 [HI (see 
Fig. la). Its norm, S IE =0.91, is close to the s-wave 
probability of 0.94 established in the deuteron. 

For closed shell nuclei, l( E {r) depends only on V eS and 
does not depend on the effective interactions determining 
their energies. Thus, the SFs for these nuclei, together 
with their ANCs, can serve in the future as a reference 
for callibrating the interaction V cS . Here, I calculate 
the overlaps ( 3 H|d), ( 4 He| 3 He) and ( 16 0| 15 N), involving 



closed shell nuclei, using M3YE. Only one Slater deter- 
minant has been used for $>a and which are divided 
by the 0s centre-of-mass motion wave function. The r osc 
is chosen to be 1.53 for 3 H and 3 He, 1.33 fm for 4 He and 
1.8 fm for 15 N and 16 to reproduce their r.m.s. radii. 
For A=3 and A=A I? E (r) are slightly smaller than the 
ab-initio overlaps from [l9|> [2(| (see Fig.l b-d) but for 
A=16 Jy (r) is slightly larger than the overlap function 
derived from the (e, e'p) knockout In both cases, S IE 
are reduced with respect to S DE (see Table I). 

For open shell nuclei, S IE also depend on occupancies 
of the single-particle orbits in the chosen model space, 
or on weights of the SU (3) and SU (4) configurations in 
the supermultiplet shell model. I generate these weights 
using phenomenological interaction from [2l| which gives 
improved spectra of Op shell nuclei. I remove the centre- 
of-mass motion explicitly and use r osc chosen as an aver- 
age of values for nuclei A and B derived in from elec- 
tron scattering. The resulting SFs S IE for ground states 
of the Op-shell nuclei, obtained with M3YE, are compared 
in Table I to S DE and to SFs available from knockout and 
those transfer reactions that use Hartree-Fock wave func- 
tions for transfer states. For all of them, S IE < S DE , 
which clearly displays the SFs reduction phenomenon. 
However, S IE > S exp for 0pi/ 2 and S IE < S exp for 
0P3/2- Agreement between S and S exp can be improved 
by tuning the V potential. In this letter, for demon- 
stration purpose only, I make the following changes to 
M3YE. All potentials in even partial waves are multi- 
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FIG. 2: (Colour online) The ratio S IE /S DE , calculated with 
RM3YE. 

plied by 1.05. This increases the SFs for d, 3 H, 3,4 He by 
10%. Then the central and spin-orbital odd components 
are multiplied by 1.7 and 2.5 respectively, which allows 
Sexp for both 12 C and 16 to be reproduced. Increasing 
odd tensor component twice reproduces the SF for 13 C. 
The SFs calculated with such a renormalised potential, 
called here RM3YE, are shown in Table I. Most SFs agree 
well with experimental data. Detailed discussion of this 
comparison will be published elsewhere. 

The ratio Rjf E = S IE /S DE , obtained with RM3YE, 
is shown in Fig. 2 as a function of AS*. The decrease 
towards large positive AS* can partially be explained by 
the presence of K-dependent function Gi(r, r') in Eq. (|3j). 
Computer calculations show that, for fixed <I>a, 3>b and 
V cS , S IE decreases with increasing n. Other effects must 
be also responsible for Rjf E (AS) behaviour but no rig- 
orous explanation to it is yet available. 

The ratio R I D E E (AS) is remarkably similar to R S (AS) 
from Q. This suggests that what really is measured in 
one nucleon removal experiments is not S DE but S IE , 
thus implying that these experiments study not occu- 
pancies of the shell model orbits but effective interac- 
tions VnCnb for occupancies fixed from other observ- 
ables, such as binding energies or nuclear spectra. Due to 
the presence of the Green function in Eq. ([3]), S IE carries 
much more information about missing model spaces than 
S DE . Therefore, it may be difficult to get correct values 
for SFs by overlapping wave functions directly even if 
they are obtained in a correlated ab-initio approach. In- 
deed, the VMC SFs for light nuclei are systematically 
larger than S IE calculated in a much simpler model with 
a reasonably chosen effective interaction, and, except for 
7 Li, the VMC SFs are in a worse agreement with exper- 
iment than those from the present work (see Table I). 

Thus, for fifty years SFs have been calculated in a pro- 
cedure of direct overlapping model wave functions that is 
sensitive only to effective interactions in truncated model 
space and does not contain important contributions from 
excluded model spaces. Calculating SFs from l[ E (r) gen- 
erated by Eq. © is a more appropriate procedure that 
allows small model spaces to be used to explain the large 
reduction of spectroscopic strength due to coupling to 



missing model spaces. Moreover, explicitely depending 
on NN matrix elements both in truncated and excluded 
spaces and having a guaranteed correct asymptotic form, 
Ii E (r) itself becomes an interface between nuclear struc- 
ture and nuclear reactions theories. Incorporating Eq. 
© into widely used shell model codes and into other 
microscopic approaches, including ab-initio ones, would 
be highly benefitial for modern nuclear physics and for 
astrophysical applications in particular. 
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